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Abstract 

We present results for the lowest order cross sections, calculated with the complete set of 
the standard model Feynman diagrams, of all possible detection channels of the associated 
production of the top quark pair and the light Higgs boson, which may be used for determi- 
nation of the top-Higgs Yukawa coupling at the future e~^e~ linear collider. We show that, 
for typical particle identification cuts, the background contributions are large. In particular, 
the QCD background contributions are much bigger than could be expected when taking 
into account a possibly low virtuality of exchanged gluons. Moreover, we include the initial 
state radiation effects and discuss the dependence of the cross sections on the Higgs boson 
and top quark masses. 
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1 Introduction 



A reaction of associated production of a top quark pair and a Higgs boson at the future 
e~^e~ linear collider 

e+e" ^ ttH (1) 

can be used to determine the top-Higgs Yukawa coupling [1]. In the Standard Model (SM), 
after taking into account decays of the top quark, t bW~^, of the antitop quark, i bW~ , 
together with the subsequent decays of the W bosons, and of the Higgs boson that, if its 
mass is less than about 140 GeV, decays preferably into a 66-quark pair, reaction ([T]) takes 
the following form 

e+e- bbbbfJihh. (2) 

If we neglect the Cabibbo-Kobayashi-Maskawa (CKM) mixing in ([2]) then /i, /2 = fe, i^n, ^r, 
u, c and f[,f2 = e~, r~, d, s, respectively. For the heavier Higgs boson, with mass uih > 
140 GeV, the decay into electroweak boson pair will dominate and the detection channels of 
([1]), will heave 10 fermions in the final state. This option is much more complicated and it 
will not be addressed in the present work. 

Reactions ([2]) receive contributions from tens thousands of Feynman diagrams already in the 
lowest order of SM. For example, reactions 



e e 



bbbbu^.e^ ^^Ufj^, 

bbbbud^^u^, 

bbbbcssc, 



(3) 
(4) 
(5) 



which correspond to the leptonic, semileptonic and hadronic decays of the W bosons, receive 
contributions, respectively, from 37868, 26 816 and 240 966 Feynman diagrams in the unitary 
gauge, with the neglect of the Yukawa couplings of the fermions lighter than the c quark. An 
overwhelming majority of the diagrams comprise background to the resonant production and 
decay of the top quark pair and Higgs boson that is mediated by 20 signal Feynman diagrams 
which contain the propagators of the top, antitop and Higgs at a time. The representative 
signal Feynman diagrams are shown in Fig. [H 




Figure 1: Representative signal Feynman diagrams of reaction ([2]) in the unitary gauge. The 
remaining diagrams are obtained by 4 permutations of the two b and two b lines. The Higgs 
boson coupling to electrons has been neglected. 
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To which extent the background contributions may affect the associated production of the top 
quark pair and Higgs boson has been discussed in [2] for a few selected reactions ([2]). A similar 
issue was also discussed in [3], where processes of the form e+e~ bbbbW^W~ bbhbl^viqq' 
accounting for the signal of associated Higgs boson and top quark pair production, as well 
as several irreducible background reactions, were studied and in |3], where pure electroweak 
(EW) contributions to the leptonic and semileptonic reactions ([2]) were computed. The off 
resonance background contributions in e^e~ Hhbud^'Vf^ have been calculated in [5]. 

In the present work, we will present results for the lowest order cross sections of all reactions 
([2]) possible in the SM calculated with the complete set of the Feynman diagrams with cuts 
on angles and energies of the final state particles that should allow for identification of the 
corresponding jets and/or separate charged leptons [6]. Such results for reaction (jl]) were 
already shown in our former work [2]. The size of the background contributions can be 
then inferred from comparison of these cross sections with the signal cross sections of the 
associated production of the top quark pair and Higgs boson. To illustrate the size of the 
QCD background we present also the cross sections calculated with the neglect of the gluon 
exchange contributions. The calculation has been performed with carlomat, a recently 
released program for automatic computation of lowest order cross sections [7]. For testing 
purposes, the lowest order cross sections calculated with the complete set of the Feynman 
diagrams have been also computed with another multipurpose Monte Carlo (MC) generator 
WHIZARD/OMEGA |S]. Moreover, we will include the initial state radiation (ISR) effects and 
discuss the dependence of the cross sections on the Higgs boson and top quark masses. 



2 Calculational details 



The results presented in Section 3 have been mostly obtained with carlomat [7], a new 
multipurpose program for automatic computation of the lowest order cross sections whose 
main goal is to provide the reliable description of multiparticle reactions. Both carlomat 
and the routines for the MC computation of cross sections it generates are written in Fortran 
90/95. 

The routine for matrix element calculation utilizes the helicity amplitude method. Poles in 
the propagators of unstable particles are regularized with constant particle widths Ta which 
are introduced through the complex mass parameters by making the substitution 

ml ^ Ml = ml- iniaTa, a = Z,W, H, t. (6) 

Substitution is made both in the s- and t-channel propagators. The calculation is per- 
formed in the complex mass scheme (CMS) [S], where the electroweak mixing parameter 
sin 6w is defined complex quantity 

sin^^v. = l-^, (7) 

with Myy and M| given by The widths of (jSD, except for Tz, whose actual value is 
rather irrelevant in the context of associated top quark pair and Higgs boson production and 
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decay, are calculated in the lowest order of SM. This, combined with the use of CMS, which 
preserves the Ward or Slavnov-Taylor identities and hence the gauge invariance, should 
minimize the unitarity violation effects. The latter are related to the fact that substitution 
(E]) introduces spurious terms of order 0{Va/ma) = 0{a) relative to the lowest order in 
places other than the resonant propagators. These higher order unitarity violating terms 
cannot be enhanced due to the exact preservation of the Ward or Slavnov-Taylor identities, 
however, the use of other widths in than the lowest order ones would have increased 
them, as we do not include any higher order corrections to the unstable particle decays. 

As in the next to leading order the number of the Feynman diagrams increases dramatically, 
calculation of the full 0{a) radiative corrections to any of reactions ([2]) is not feasible either at 
the moment or even in the foreseeable future. However, precision of theoretical predictions for 
such reactions can be improved by the inclusion of the ISR effects in the leading logarithmic 
(LL) approximation in the structure function approach. The most recent version of carlomat 
allows to compute the ISR corrected cross section in LL approximation da^^{pi,p2) according 
to the following formula 



where pi (^2) is the four momentum of a positron (electron), Xi (X2) is the fraction of 
the initial momentum of the positron (electron) that remains after emission of a collinear 
photon and d(TBorn(a^iPi, a^2P2) is the lowest order cross section calculated at the reduced 
four momenta of the positron and electron. The structure function T^^ [x, Q"^) is given by 
Eq. (67) of |10j, with 'BETA' choice for non- leading terms. The splitting scale Q^, which is 
not fixed in the LL approximation is selected to be equal s = (pi +^2)^- 

carlomat generates also dedicated phase space parametrizations which take into account 
mappings of peaks in the matrix element caused by propagators of massive unstable parti- 
cles, of a photon, or a gluon in each Feynman diagram. This means that a number of different 
phase space parametrizations generated is equal to a number of the Feynman diagrams. The 
phase space parametrizations are automatically implemented into a multichannel MC inte- 
gration routine that performs integration over a 20-dimensional phase space of reactions 
The integration is performed in several iterations, with the weights of different kinematical 
channels determined anew after each iteration. In this way, the kinematical channels with 
small weights are effectively not used in the subsequent iterations. 



3 Results 

We use the following set of initial physical parameters: the Fermi coupling, fine structure 
constant in the Thomson limit and strong coupling 

= 1.16639 X 10"^ GeV-^ ao = 1/137.0359991, a,(mz) = 0.1176, (9) 

the W- and Z-boson masses 



dcTBorn+ISR 




JO Jo 



mw = 80.419 GeV, 



mz = 91.1882 GeV, 



(10) 
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the top quark mass and the heavy external fermion masses 



mt = 174.3 GeV, = 4.8 GeV, m, = 1.3 GeV, = 1.77699 GeV. (11) 

The same masses are used both in the matrix elements and phase space parametrizations. 
We would like to stress that, in the CMS, the complex top quark mass calculated according 
to Eq. enters both the top quark propagator, including the numerator, as well as the 
top-Higgs Yukawa coupling. In order to speed up the computation fermion masses smaller 
than those of the c quark and r lepton have been neglected. 

The value of the Higgs boson mass is assumed at mn = 130 GeV. To avoid unitarity violation 
the widths of t-quark, W- and Higgs bosons are calculated to the lowest order of SM resulting 
in the following values: 

Tt = 1.53088 GeV, Tw = 2.04764 GeV, Th = 8.0540 MeV. (12) 

The Z boson width, whose actual value is not relevant in the calculation, is fixed at its 
experimental value Tz = 2.4952 GeV. 

We identify jets with their original partons and define the following basic cuts which should 
allow to detect events with separate jets and/or isolated charged leptons: 

• cuts on an angle between a quark and a beam, an angle between two quarks and on a 
quark energy: 

5° < 0(g,beam) < 175°, e{q,q')> 10°, > 15 GeV, (13) 

• cuts on angles between a charged lepton and a beam, a charged lepton and a quark 
and on energy of the charged lepton, / = /x, r, in the semileptonic detection channels 

5° < ^(/,beam) < 175°, ^(/,g)>10°, > 15 GeV, (14) 

• a cut on an angle between the two charged leptons in the leptonic detection channels 
of® 

^(/,/')>10°, (15) 

• a cut on the missing transverse energy in the hadronic and semileptonic detection 
channels of ([2]) 

> 15 GeV. (16) 

The results for the lowest order cross sections of the different channels of ([2]), which corre- 
spond to different decay modes of the and W~ bosons resulting from t and i decays, 
are shown in Tables [lH3] for a/s = 500 GeV, 800 GeV, 1 TeV and 2 TeV, the first two centre 
of mass energies being characteristic for the International Linear Collider (ILC) project [H] 
and the second two for the Compact Linear Collider (CLIC) design [12]. We present cross 
sections of specific reactions, calculated with the complete set of the lowest order Feynmamn 
diagrams: cx^'^'^-, computed with WHIZARD/OMEGA and CTaii, computed with carlomat. Both 
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(T^ and (Tall have been computed with the same cuts and initial input parameters. The 
comparison shows agreement within 1-3 standard deviations of the MC integration that is 
indicated in parentheses. In order to illustrate the size of the QCD background for each 
reaction we present also the cross sections cr^o qcd with the neglect of the gluon exchange 
contributions, computed with carlomat. Comparison of cXaii and (TnoQCD shows that the 
QCD contributions are sizeable, amazingly enough also for the leptonic channels of Tabled! 
They are particularly large at a/s = 500 GeV, where the signal is reduced because of the 
limited phase space volume available for reaction ([1]). However, the QCD background is 
sizeable also for higher centre of mass energies. This is somewhat surprising if one takes into 
account a possibly low virtuality of the exchanged gluons. Let us note that, despite quite 
different numbers of the Feynman diagrams contributing, the cross sections in Tables [TH3] 
are practically flavor independent within each of the detection channels. This means that, 
due to the cuts enforced, there is practically no dependence on m^, rric and m,- in the gluon 
and/or photon exchange background contributions. 

The full size of the background contributions can be inferred from comparison of the cross 
sections (Xaii of Tables [TH3] with the corresponding signal cross sections presented in Table HI 
The background contributions are large, but they can be efficiently reduced by imposing 
cuts on invariant masses allowing for reconstruction of the top, antitop and Higgs boson. 
In particular, a narrow cut, of the order of 1 GeV, on the invariant mass of the two 6-jets 
that reconstruct the Higgs boson mass reduces the background very efficiently. This issue 
has been already extensively discussed in [2j so we will not address it here. The signal cross 
sections are also the same for all the channels within a given type: leptonic, semileptonic 
and hadronic one, as we have neglected the light fermion masses. However, even if the light 
fermion masses had been kept non zero, the flavor dependence of the cross sections would 
have been negligibly small, as amplitudes of the signal diagrams of the associated production 
of the top quark pair and Higgs boson are almost completely independent of the external 
fermion masses which always couple to the propagators of the heavy particles. 

Let us note that, in spite of somewhat different cuts imposed, the cross sections of the 
semileptonic channels in Tables [THU are roughly a factor 3 bigger than those of the leptonic 
channels. A similar relation holds between the cross sections of the hadronic and semileptonic 
channels. 

It is interesting to look at the dependence of the cross sections on the Higgs boson and 
top quark masses. In order to illustrate this we show in Table [5] the cross sections Uaii, 
c^no QCD and (Tsignai of rcactiou (jl]) with cuts given by f|T3l) . f|T^ and f|T6|) that have been 
calculated with 2 different hypothetical values of the Higgs boson mass: mn = 115 GeV 
and rriH = 130 GeV, and 3 values of the top quark mass: rrit = 171 GeV, rrit = 174.3 GeV 
and rrit = 176 GeV. The top quark mass has been chosen within roughly 2a range of the 
combined Tevatron result rrit = 173.1 ±1.3 GeV [13]. The change of particle masses affects 
the corresponding lowest order SM widths of (fT2l) which become 

Tt{171 GeV) = 1.43152 GeV, Tt{l76 GeV) = 1.58351 GeV, 

1^(115 GeV) = 6.0223 MeV. (17) 

The cross sections for rriH = 115 GeV are bigger than those for rriH = 130 GeV. The relative 
difference is largest at y/s = 500 GeV, where it amounts to a factor 3-4 for the signal cross 



6 





e+e — > hhhhve^'^ 11 


e+e —>■ bbbbVf.e^T 


No. of graphs: 


37 868 33 500 


42 428 37 708 


[GeV] 


_Whiz. 

<7all '^all CTnoQCD 


_Whiz. 

O^all CTall CTnoQCD 


500 
800 
1000 
2000 


8.71(15) 8.26(6) 2.889(7) 
36.2(1) 35.6(2) 24.41(4) 
M 1 (W 34 Hfo^ 22 50f4^ 
18.6(2) 18.4(2) 11.15(7) 


8.47(5) 8.28(6) 2.882(7) 
36.1(1) 36.0(1) 24.42(4) 
MliW 34 3fl^ 22 57f5') 
18.2(2) 18.3(1) 11.06(7) 








No. of graphs: 


131 648 115 200 


46 890 40850 


[GeV] 


_Whiz. 

C^all '^all CTnoQCD 


_Whiz. 

O^all CTall CTnoQCD 


500 
800 
1 000 
2000 


8.43(4) 8.43(6) 2.895(8) 
36.1(1) 36.1(1) 24.57(4) 
34 2fl^ 34 3fl^ 22 55f4^ 
19.3(2) 20.0(2) 12.28(14) 


8.43(1) 8.49(7) 2.890(8) 
35.8(1) 35.9(1) 24.48(4) 
34 0fl^ 22 54f4') 
17.6(1) 17.7(1) 10.52(3) 




e'^e~ —>■ bbbburT^T^Ur 


e+e^ ^ bbbbi'fj^fi'^T'Ur 


No. of graphs: 


66 658 58 666 


21 214 18 854 


[GeV] 


_Whiz. 

O^all CTall CTnoQCD 


_Whiz. 

C^all CTall CTnoQCD 


500 
800 
1000 
2000 


8.44(1) 8.43(5) 2.907(6) 
35.8(1) 35.8(1) 24.49(4) 
34.0(2) 34.1(1) 22.48(4) 
17.5(1) 17.7(1) 10.50(3) 


8.40(3) 8.46(3) 2.885(6) 
36.2(2) 35.9(1) 24.49(4) 
33.8(1) 34.2(1) 22.42(4) 
17.3(1) 17.7(1) 10.56(4) 



Table 1: Cross sections in ab of leptonic detection channels of ([2]) calculated with 
WHIZARD/OMEGA (first column) and carlomat (second column) with cuts given by f |T3i) . ( |T^ 
and (fT5|) . The numbers in parenthesis show the MC uncertainty of the last decimal. 

section, a factor 2 for pure EW contributions and a factor 1.3-1.5 for the full cross section. 
The change of a top quark mass in Table [5] is smaller than that of the Higgs mass, however, 
it has proportionally the same effect on the cross section at ^/s = 500 GeV. This can be 
easily traced back to the phase space volume limitations that are much stronger for the 
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e+e —>■ bbbbude z/g 


e^e —>■ bbbbcse v^. 


No. of graphs: 


53 632 39 224 


59 144 43 816 


[GeV] 


_Whiz. 

<7all all (7 no QCD 


_Whiz. 

<7all '^all CTnoQCD 


500 
800 
1000 
2000 


26,8(1) 26.9(1) 7.91(2) 
99,9(3) 100.6(3) 67.3(1) 
94 ()(^) 95 61 9(1) 
49.6(7) 49.1(3) 29.4(1) 


26.8(1) 26.9(1) 7.93(2) 
99.7(5) 100.5(3) 67.1(1) 
Q2 4^4^ 94 61 '^iW 
47.8(4) 48.2(3) 29.8(2) 




e+e^ hbbbud^^v^ 


e+e^ bbbbcsfi^u^ 


No. of graphs: 


26 816 19 612 


29 572 21 908 


[GeV] 


CTall '^all (^noQCD 


^Whiz. ^ ^ 

C^all ^ all no QCD 


500 
800 
1000 

-L W WW 

2000 


26.6(1) 26.9(1) 7.91(2) 
98.6(3) 100.6(2) 67.0(1) 
93 3('3') 94 Ai?,\ 61 5f1') 

ty^.^lO/ WX.tJlXI 

46.7(2) 47.6(2) 28.3(1) 


26.5(1) 26.6(1) 7.89(2) 
98.4(3) 100.0(3) 67.3(1) 
92 4^5^ 93 61 5f1^ 
47.0(1) 46.9(2) 28.2(1) 




e^e~ — *■ bbbbudr^ Vr 


e^e' ^ bbbbcsT^ u-r 


No. of graphs: 


29 480 21 908 


32 332 24 300 


[GeV] 


_Whiz. 

C^all ^ all CT no QCD 


_Whiz. 

C^all all Cr no qcd 


500 
800 

1000 
2000 


26.8(1) 27.0(1) 7.89(2) 
99.4(6) 100.7(2) 67.3(1) 
93.5(3) 94.1(2) 61.6(1) 
47.4(3) 47.3(2) 28.3(1) 


26.6(1) 26.6(1) 7.90(2) 
100.3(2) 100.9(3) 67.3(1) 
92.8(2) 94.0(3) 61.5(1) 
47.2(1) 47.0(2) 28.3(1) 



Table 2: Cross sections in ab of semileptonic detection channels of ([2]) calculated with 
WHIZARD/OMEGA (first column) and carlomat (second column) with cuts given by f |T3i) . (HM 
and (fT6|) . The numbers in parenthesis show the MC uncertainty of the last decimal. 

heavier Higgs boson and top quark. At higher energies, where the relative differences in 
cross sections are smaller, the mass dependence becomes more involved as the phase space 
effect interferes with the shift of the maximum of the cross section dependence on ^/s that is 
caused by changes in rriH and/or rrit. For example, the signal cross section at y/s = 800 GeV 
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—>■ hhhhuddu 




e bbbbcssc 


No. of graphs: 




185 074 


81 250 




240 966 


114 190 


Vs [GeV] 


_ vv niz. 


0"all 


0" no QCD 


— . VV illZ. 

^all 


0"all 


O'no QCD 


500 


92.8(4) 


94.1(3) 


24.10(7) 


93.5(6) 


93.1(2) 


24.16(5) 


800 


318(2) 


314(1) 


205.7(3) 


308(2) 


314(1) 


205.7(3) 


1000 


284(2) 


291(1) 


187.4(3) 


286(2) 


291(1) 


187.6(3) 


2000 


137.1(5) 


139.2(7) 


83.1(3) 


137(1) 


139.2(5) 


83.0(3) 



e^e — s> bbbbudsc 



No. of graphs: 




39 342 


25 246 


[GeV] 


^Whiz. 
"^all 


0"all 


O'no QCD 


500 
800 
1000 
2000 


92.2(3) 
311(2) 
288(2) 
139(1) 


93.2(2) 
314(2) 
290(1) 
139(1) 


24.10(5) 

205.6(3) 
187.2(3) 
82.8(2) 



Table 3: Cross sections in ab of hadronic detection channels of ([2]) calculated with 
WHIZARD/OMEGA (first column) and carlomat (second column) with cuts given by ( fT3l) and 
( 1161) . The numbers in parenthesis show the MC uncertainty of the last decimal. 



v/i [GeV] 


C^leptonic 


C"semileptonic 


^hadronic 


500 


1.12(1) 


3.095(3) 


9.40(1) 


800 


16.86(4) 


46.27(2) 


141.4(1) 


1000 


14.67(4) 


40.18(2) 


122.5(1) 


2000 


5.67(7) 


15.14(3) 


44.58(3) 



Table 4: Signal cross sections in ab of different detection channels of ([2]). The cuts are given 
by ([13]), ([H]) and ([IS]) for the leptonic, by ([13]), ([H]) and ([16]) for the semileptonic and by ([13]) 
and f[T6]) for the hadronic channels. The numbers in parenthesis show the MC uncertainty 
of the last decimal. 

for tuh = 115 GeV is about 50% bigger than that for mu = 130 GeV, while the cross section 
for mi = 171 GeV is smaler smaler than that for = 176 GeV. 
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= 115 GeV 




r = 130 GeV 


[GeV] 


mt [GeV] 


0"all 


0" no QCD 


(^signal 


0"all 


<7no QCD 


'^signal 


500 


171 


39.8(1) 


20.21(5) 


14.61(1) 


29.9(1) 


10.29(3) 


4.670(4) 




174.3 


35.3(1) 


16.41(4) 


11.69(1) 


26.9(1) 


7.91(2) 


3.095(3) 




176 


33.3(1) 


14.52(4) 


10.19(1) 


25.4(1) 


6.75(2) 


2.365(2) 


800 


171 


118.2(2) 


84.8(1) 


63.17(3) 


99.2(2) 


66.0(1) 


44.62(2) 




174.3 


119.9(3) 


86.8(1) 


65.61(3) 


100.6(2) 


67.0(1) 


46.27(2) 




176 


120.3(3) 


88.0(1) 


66.83(3) 


100.4(2) 


67.8(1) 


47.08(2) 


iUUU 


III 


iOo.ol^o j 


VO O / 1 \ 


Oi. I o\z, j 


92.1(2) 


o9.o(^i ) 


oo OA 




174.3 


108.6(3) 


75.7(1) 


54.22(3) 


94.4(3) 


61.5(1) 


40.18(2) 




176 


110.2(3) 


77.7(1) 


55.51(3) 


94 


62.3(1) 


41.07(2) 


2000 


171 


49.9(2) 


31.37(9) 


18.17(5) 


46.2(2) 


27.31(7) 


14.34(4) 




174.3 


51.5(2) 


32.64(9) 


19.20(5) 


47.6(2) 


28.35(8) 


15.14(3) 




176 


51.8(2) 


33.21(9) 


19.78(5) 


47.5(2) 


28.75(9) 


15.60(3) 



Table 5: Cross sections in ab of (jlj) calculated with carlomat with different values of the 
top quark mass mt and the Higgs boson mass ttih with cuts given by (fT3|) . (fl^ and (fT6|) . 
The numbers in parenthesis show the MC uncertainty of the last decimal. 

In order to illustrate the potential effect that the unitarity violating terms discussed in 
Section 2 may have on the presented results we show in Table [6] the lowest order cross 
sections of (jlj) calculated with carlomat for two different values of the Higgs boson mass 
niH and mt = 174.3 GeV in the CMS and fixed width scheme (FWS). The latter differs from 
the CMS in the definition of the EW mixing parameter that instead of by Eq. is defined 
by 

sin2^VF = l-^, (18) 

with mw and mz being physical masses of the EW bosons. As the cross sections agree 
within one standard deviation of the MC integration, we see that the unitarity violation 
terms are not enhanced in the FWS at the presented values of ^/s, even though the Ward 
or Slavnov-Taylor identites are not stisfied in this scheme. 

The impact of the ISR is illustrated in Table [3, where we show the cross sections (Xaii+isR, 
c"no QCD+isR and cTsignai+isR ^01 reaction (jl]) including ISR according to Eq. ([8]) with the 
splitting scale = s. They have been calculated with the Higgs boson masses mn = 
115 GeV and mn = 130 GeV, and the top quark mass mt = 174.3 GeV. These cross sections 
should be compared with the Born cross sections at the corresponding centre of mass energy 
and mt = 174.3 GeV in Table [51 We see that, as could be expected from Eq. ([8]), ISR 
decreases substantially the signal and pure EW cross sections at y/s = 500 GeV, as due to 
the photon emission the actual centre of mass energy of e^e~ scattering may fall down below 
the tiH threshold, where the signal cross section cisignai becomes very small. Note that the 
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niH = llh GeV 


niH = 130 GeV 


[GeV] 




(Tail cr^n 


500 
800 
1000 
2000 


35.3(1) 35.4(1) 
119.9(3) 120.3(3) 
108.6(3) 108.7(2) 
51.5(2) 51.3(2) 


26.9(1) 26.8(1) 
100.6(2) 100.8(3) 
94.4(3) 94.5(2) 
47.6(2) 47.3(2) 



Table 6: Cross sections in ab of (jl]) in the CMS and FWS calculated with carlomat for two 
different values of the Higgs boson mass rnn and = 174.3 GeV. Cuts are given by (IT^ . 
(HM and (fTB]l . The numbers in parenthesis show the MC uncertainty of the last decimal. 

effect is smaller for itlh = 115 GeV due to more phase space volume available. The decrease 
of the full cross section at this energy is relatively smaller, as it is dominated by the QCD 
background contributions. At y/s = 800 GeV and y/s = 1 TeV, the cross section reduction 
is much smaller and, at y/s = 2 TeV, there is even some growth in the ISR corrected cross 
sections with respect to the lowest order cross sections of Table O 





= 115 GeV 


mH = ISO GeV 


[GeV] 


(^all+ISR 


0"no QCD+ISR 


0"signal+ISR 


0"all+ISR 


O"no QCD+ISR 


0"signal+ISR 


500 


29.2(1) 


12.24(4) 


8.48(2) 


22.6(1) 


5.89(2) 


2.114(6) 


800 


114.1(2) 


82.2(1) 


62.34(6) 


94.3(2) 


62.9(1) 


43.39(4) 


1000 


107.5(3) 


75.3(1) 


54.29(5) 


91.9(2) 


60.5(1) 


39.73(4) 


2000 


54.7(2) 


35.2(1) 


21.17(5) 


50.0(2) 


30.2(1) 


16.50(3) 



Table 7: Cross sections in ab of (jlj) including ISR in the LL approximation calculated with 
carlomat for two different values of the Higgs boson mass rriH and rrit = 174.3 GeV. Cuts 
are given by (fT3|) . (HM and (fT6l) . The numbers in parenthesis show the MC uncertainty of 
the last decimal. 



4 Summary 

We have presented results for the lowest order cross sections, calculated with the complete 
set of the standard model Feynman diagrams, of all reactions ([2]) relevant for the associated 
production of the top quark pair and Higgs boson that can be used for determination of the 
top-Higgs Yukawa coupling at the e~^e~ linear collider. A comparison with the corresponding 
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signal cross sections of ([T]) has shown that the background contributions are large for typical 
particle identification cuts. In particular, the QCD background contributions are much 
bigger than it could have been expected taking into account a possibly low virtuality of 
exchanged gluons. As we have shown elsewhere [2] for a few representative channels of ([2]) 
the background can be efficiently reduced by imposing invariant mass cuts allowing for the 
top and antitop quark, and Higgs boson identification. We have shown that the unitarity 
violating terms discussed in Section 2 have practically no effect on the lowest order cross 
sections of (jl]) at the considered range of y/s by comparing the cross sections calculated in 
the CMS and FWS, where the latter violates gauge invariance. 

Moreover, we have included the ISR and illustrated its effects concentrating on a semileptonic 
reaction (jlj) for which we have have also discussed the dependence of the cross sections on 
the Higgs boson and top quark masses. Taking into account cross sections of the different 
tiH detection channels presented in Tables [TH7] the best place to measure the top-Higgs 
Yukawa coupling seem to be a linear collider operating at the centre of mass energy of about 
800 GeV. 
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